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We extend our previous studies of the neutrino fluxes expected from neutralino 
LSP annihilations inside the Sun to include variants of the minimal supersymmetric 
extension of the Standard Model (MSSM) with squark, slepton and gaugino masses 
constrained to be universal at the GUT scale, but allowing one or two non-universal 
super symmetry-breaking parameters contributing to the Higgs masses (NUHM1,2). 
As in the constrained MSSM (CMSSM) with universal Higgs masses, there are large 
regions of the NUHM parameter space where the LSP density inside the Sun is not 
in equilibrium, so that the annihilation rate may be far below the capture rate, 
and there are also large regions where the capture rate is not dominated by spin- 
dependent LSP-proton scattering. The spectra possible in the NUHM are qualita- 
tively similar to those in the CMSSM. We calculate neutrino-induced muon fluxes 
above a threshold energy of 10 GeV, appropriate for the IceCube/DeepCore detector, 
for points where the NUHM yields the correct cosmological relic density for repre- 
sentative choices of the NUHM parameters. We flnd that the IceCube/DeepCore 
detector can probe regions of the NUHM parameter space in addition to analogues 
of the focus-point strip and the tip of the coannihilation strip familiar from the 
CMSSM. These include regions with enhanced Higgsino-gaugino mixing in the LSP 
composition, that occurs where neutralino mass eigenstates cross over. On the other 
hand, rapid-annihilation funnel regions in general yield neutrino fluxes that are un- 
observably small. 
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I. INTRODUCTION 



One of the most actively pursued strategies for detecting supersymmetric dark matter 
particles (XSPs) is the search for signatures of the annihilations of LSPs inside the Sun 
or Earth [J, The principle for detection is to search for the passage through a large 
detector in ice or water of muons produced by the interactions of energetic neutrinos released 
in the LSP dark matter annihilation process. There have been extensive studies of this 
potential experimental signature in many variants of the MSSM, and experiments such as 
IceCube/DeepCore 0-0] are starting to chip away at the MSSM parameter space. 

We recently re-analyzed 0] this potential signature in the framework of the MSSM with 
all the supersymmetry-breaking spin-1/2 and -0 mass parameters {mi/2,mQ) constrained to 
be universal at the GUT scale (the CMSSM) [8i], imposing the requirement that the LSP 
should provide the density of dark matter inferred from WMAP and other experiments. 
In the CMSSM, where the supersymmetry breaking trilinear mass parameters, Aq are also 
taken to be universal at the GUT scale, the resulting relic density is found to lie in the 
WMAP range only along relatively narrow strips in the (^«i/2, '"^o) plane for fixed tan /3 and 



Aq [10|, lUj. These correspond to the coannihilation strip, where the mass of the lightest 
neutralino is close to the mass of the lightest charged slepton (usually the mostly right- 
handed stau); the heavy Higgs funnel, found at large tan/3 and large {mi/2,mo), where the 
neutralino mass is close to half the heavy Higgs mass and rapid annihilations of neutralinos 
are mediated by the s-channel exchange of heavy Higgs scalars and pseudoscalars; and the 
focus-point region which is typically found at very large values of mg when the /x parameter 
(an output of the minimization of the Higgs potential in the CMSSM) is driven to small 
values and the neutralino picks up a more significant Higgsino component. 

In our previous work |/7,], we found that the LSP capture rate was not in general dominated 
by scattering on protons inside the Sun via spin-dependent couplings, but that an important 
role was often played by spin-independent scattering on heavier nuclides. We also found that, 
in many regions of the CMSSM parameter space, LSP capture and annihilation would not 
be in equilibrium, and that the annihilation rate would be correspondingly reduced. We 
also analyzed the uncertainties in the magnitude of the potential muon-neutrino signal due 
to uncertainties in the composition of the Sun and in the scattering matrix elements. We 
found that the CMSSM might be detectable in IceCube/DeepCore along (some part of) the 
WMAP strip in the focus-point region of the (mi/^, mo) plane, and near the low-mi/2 tip of 
the WMAP strip in the coannihilation region 0, 12|. 

In this paper we extend these previous studies to models with one or two degrees of 
non- universality in the soft supersymmetry-breaking contributions to the Higgs doublets, 
the NUHMl [3, [l| and NUHM2 . One of our primary objectives is to understand 



the circumstances under which such relatively high neutrino fluxes may be attained in these 
models and, conversely, whether the relatively low rates usually found in the CMSSM are 
specific to that model. More generally, we seek to lay a basis for systematic comparisons of 
the physics capabilities of different detection strategies in (relatively) simple variants of the 
MSSM. 



Supersymmetric dark matter searches [6|, Il7|, Il8| are complementary to searches for super- 
symmetry at accelerators [l^, 20 1. At the moment, the latter are sensitive primarily to the 
spin-1/2 and -0 mass parameters mi/2 and mo, and are less sensitive to the non-universality 
parameters that appear in the NUHMl, 2. For example, global likelihood fits j2l| currently 
yield similar 68 and 95% confidence- level preferred regions in the (mi/2,mo) planes of the 
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CMSSM and NUHMl. It is therefore particularly interesting to know whether direct and 
indirect dark matter searches offer ways to differentiate between these models. 

The NUHMl, 2 offer additional mechanisms to bring the relic LSP density into the WMAP 
range, in addition to the coannihilation, rapid-annihilation and focus-point possibilities men- 
tioned above in the context of the CMSSM. For example, there are distinctive regions of 
NUHM parameter space where Higgsino-gaugino mixing in the LSP is enhanced by level- 
crossing in the neutralino mass matrix, bringing the relic density into the WMAP range. 
Alternatively, the LSPs may annihilate rapidly through direct-channel heavy Higgs if, A 
poles even if mi/2,'^o and tan/3 are relatively small. It was shown in [2l[ that the col- 
lider prospects for sparticle detection are rather different in the NUHMl low-mass rapid- 
annihilation region than they are in the CMSSM, whereas the favoured rates for direct LSP 
detection via scattering on nuclei were broadly similar in the CMSSM and the NUHMl 
(though the uncertainties were greater in the latter case). Therefore, it is interesting to 
study the prospects in this region for LSP detection via the energetic neutrinos produced 
by annihilation inside the Sun. 

The layout of this paper is as follows. In Section [IT] we recall briefly some general features 
of the NUHMl, 2 and discuss other inputs into the rate calculations. Then, in Section UTTl we 
explore the solar annihilation rates in some generic slices through the NUHMl parameter 
space, finding that they are enhanced in regions with relatively large Higgsino components 
in the LSP as may occur for specific relations between /i and mi/2 where there is level 
crossing. In Section IIVI we extend our analysis to slices through the NUHM2 parameter 
space. Section IVl summarizes our conclusions. 

We find that whereas the neutrino flux may be observable in regions with enhanced 
Higgsino-gaugino mixing, analogously to the enhancement along the focus-point WMAP 
strip in the CMSSM, the flux is generically unobservably low in the rapid-annihilation fun- 
nels. This suggests that the observation of a high-energy solar neutrino flux in the Ice- 
Cube/DeepCore experiment is a potential diagnostic for large mixing and level crossing in 
the neutralino mass matrix, and specifically of the relation between fi and mi/2, which is a 
potential tool for identifying non-universal Higgs mass parameters. 



II. PREAMBLE 

A. The NUHMl and NUHM2 Parameter Spaces 

In the CMSSM, the free parameters are the supposedly universal supersymmetry-breaking 
parameters mi/2, ruQ and Aq, as well as tan (3. The Higgs mixing superpotential parameter n 
and the bilinear supersymmetry-breaking parameter Bq, and hence the pseudoscalar Higgs 
mass rriA, are then determined using the electroweak vacuum conditions, with a sign ambi- 
guity in /i. The sign of the discrepancy between the experimental value of — 2 and the 
value calculated within the Standard Model suggests that > 0, and CMSSM analyses are 
often presented in {mi/2,mQ) planes for yU > and fixed values of tan/3 and Aq. Although 
much of our NUHM analysis is for > 0, we also consider the possibility of a negative sign. 



The value of Aq is notoriously unconstrained, see, e.g., [22|, and for definiteness we set it to 
zero in what follows. 

In the NUHMl, the soft supersymmetry-breaking contributions to the masses of the two 
MSSM Higgs doublets, mi and m2, are assumed to be equal, but are allowed to differ from 
rriQ. The extra degree of freedom may be used to treat either or as a free parameter. 
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On the other hand, in the NUHM2, the two soft supersymmetry-breaking contributions to 
the masses of the MSSM Higgs doublets are allowed to vary independently, and the two extra 
degrees of freedom may be used to treat both and as free parameters. Analyses of the 
NUHMl and NUHM2 parameter spaces are often presented in planes spanned by pairs of 
the quantities mi/2,mo,mA and fi for some fixed values of tan/3 and the other parameters, 
and we present some examples below which are selected from the analysis in Since 
a frequentist likelihood analysis favours relatively small values of mi/2 and mo in both 
the CMSSM and NUHMl Hll, we concentrate here on NUHMl and NUHM2 planes with 
relatively low values of either mi/2 and mo- We note, however, that a complete likelihood 
analysis of the NUHM2 is yet to be performed. 



B. Spin-Dependent and -Independent Scattering Rates 

When calculating the LSP annihilation rates inside the Sun, the key particle physics 
inputs — apart from the choice of supersymmetric model — are the matrix elements for dark 
matter scattering on the nuclides inside the Sun. It is often assumed that LSP capture in 
the Sun is dominated by spin-dependent scattering on hydrogen but, as discussed in 
spin-independent scattering on heavier nuclei actually dominates in generic regions of the 
CMSSM parameter space. Figure [T] displays contours of the ratio of the solar dark matter 
annihilation rate calculated using only spin-dependent scattering to the total annihilation 
rate including also spin-independent scattering in (left) the CMSSM (mi/2,'^o) plane for 
tan/? = 10, Aq = and /i > 0, and (right) the NUHMl {mi/2,mo) plane for tan/3 = 10, 
Aq = and /i = 500 GeV. Regions excluded because there is no consistent electroweak 
symmetry breaking (EWSB) have dark pink shading, those with a charged LSP have brown 
shading, and those in conflict with 6 — )• 57 measurements have green shading. Regions to 
the left of the black dashed (red dash-dotted) line are inconsistent with the absence at LEP of 



charginos (a Higgs boson) |24j-|26|. We recall that there is a theoretical uncertainty ~ 1.5 GeV 
in the calculation of rrih in the CMSSM, which induces an uncertainty ~ 50 GeV in values of 
mi/2 along the red dash-dotted line. The recent LHC direct exclusion regions in the CMSSM 



nT'i/2, iTT'o) plane [19|, |20[ are somewhat weaker than the indirect LEP Higgs constraint shown 



here. The pale pink band is favoured by the BNL measurement of — 2 |27H29|. at the 
±1(2) — a level along the dashed (solid) lines, but we do not impose this as a constraint 
on our analysis. In the CMSSM case (left panel), we see that spin-dependent scattering is 
dominant only at small mi/2 and large mo. In the NUHMl case (right panel), we see that 
spin-dependent scattering is subdominant at small mi/2, becoming more important along 
the WMAP-compatible strip with mi/2 ~ 1000 GeV. However, spin-independent scattering 
is important even here, and in the following analysis all calculations of the annihilation rates 
include the contributions from both spin-dependent and -independent scattering. 

The uncertainties in the spin-independent scattering matrix element include those in the 
ratios of the light quark masses, the octet contribution, (Tq, to the pion-nucleon cr term, Sjrjv, 
and the value of Sj^jv itself [30H34j] . The largest uncertainty is due to S^jv, for which we used 
the value 64 MeV as our default in whilst also exploring the implications of other values. 
The second-largest uncertainty is due to that in ctq, for which we used the value 36 MeV as our 
default. We assume the same default values in this analysis. As discussed in more detail in 
Refs. 0,|33[, different measurements for E^riv and do lead to variations in the spin- independent 
scattering cross-section by a factor of ~2-3 and the choice of experimental values for these 
two parameters significantly impacts the spin-independent scattering contribution to capture 
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mi/2 (GeV) mi/2 (GeV) 



FIG. 1: Sample {1111/2, nT-o) planes for tan/3 = 10 and Aq = in (left) the CMSSM for /i > 
and (right) the NUHMl for jj, = 500 GeV, showing regions excluded because there is no consistent 
electroweak vacuum (dark pink shading), or because there is charged dark matter (brown shading), 
or because of a conflict with 6 — ?> 37 measurements (green shading). Only regions to the right of the 
black dashed (red dash-dotted) line are consistent with the absence at LEP of charginos (a Higgs 
boson). The turquoise strips are favoured by the determination of the cold dark matter density by 
WMAP and other experiments \^], and the light pink band is favoured by the BNL measurement 
of 9 II — 2. We also show contours of the ratio of the solar dark matter annihilation rate calculated 
using only spin- dependent scattering to the total annihilation rate including also spin-independent 
scattering. 

in the Sun. If St^jv were smaller, the total annihilation rate would decrease, and the ratio 
of the solar dark matter annihilation rate calculated using only spin-dependent scattering 
to the total annihilation rate including also spin-independent scattering shown in Figure [1] 
would increase. By comparison, the uncertainties in the light quark mass ratios are much 
less significant. 

The principal uncertainty in the spin- dependent matrix element is the contribution of 
strange quarks to the nucleon spin, Ag, with the uncertainties due to qa and the SU(3)- 
octet nucleon matrix elements being significantly smaller. The range —0.06 > As > —0.12 
was considered in and was found to induce an uncertainty in the annihilation rate that 
was ~ 10%, considerably smaller than the others considered. In the following we take the 
central value Ag = —0.09 [ssl], which was adopted in 0] as the default value. 



C. Capture/ Annihilation Rates and Neutrino/Muon Fluxes 

We compared in [3] the rates estimated in various alternative solar models. We found 
that there was at most a 4% difference in the annihilation rates between the two models 
of Serenelli et al. [s^ (AGSS09 and AGSSOQph) based upon recent abundance estimates 



37[, and assumed the AGSS09 model as our default. We do the same here, performing a 
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mi/2 (GeV) mi/2 (GeV) 

FIG. 2: The same CMSSM plane (left) and NUHMl plane (right), displaying also contours of the 
ratio of solar dark matter annihilation and capture rates. Equilibrium corresponds to a ratio of 
unity, which is approached for small mi/2 o.'^d large mo in the CMSSM, and near a vertical strip 
with mi/2 ~ 1000 GeV in the NUHMl example. 

full numerical integration over the radial profile of the Sun when determining the capture 
rates ^. 

In modelling the dark matter halo, we assume a non-rotating isothermal sphere with an 
rms speed of 270 km/s, a disk rotation speed of 220 km/s, and a local dark matter density of 
0.3 GeV/cm^. If the calculated neutralino relic density is below the WMAP observed relic 
density, we assume that only a fraction of the local dark matter density, equal to the ratio 
of the neutralino and WMAP dark matter relic densities, is attributable to neutralinos. 
As in 0], we do not address other halo models in this paper, but note that our results 
would scale linearly with the local dark matter density as equilibrium between capture and 
annihilation is approached. 

However, as discussed extensively in 0], equilibrium is not in general reached in the 
CMSSM, and an example is shown in the left panel of Figure [21 namely the (^^1/2, ^0) plane 
for tan /3 = 10, = and /i > 0. The dark blue contours are for different ratios of the 
annihilation and capture rates, and we see that equilibrium is reached only for small mi/2 
and large ttiq. The right panel of Figure [2] shows the corresponding contours in the NUHMl 
(mi/2, itiq) plane for tan (3 = 10, Aq = and fi = 500 GeV, where we see that equilibrium is 
approached only near a vertical strip with mi/2 ~ 1000 GeV. In the following we calculate 
annihilation rates without relying on the assumption of equilibrium. 

The neutralino annihilations produce high-energy neutrinos which, through interactions 
with matter, will induce muons in or around a detector such as IceCube/DeepCore 
The IceCube detector has outfitted ~1 km^ of ice at the South Pole with optical sensors 



Using the simpler Gould approximation [SSl would have yielded rates differing from the exact results by 
at most 6%, as was shown in |7|]. 
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to observe the Cerenkov light produced by the passage of muons through the ice. The 
large volume allows for sensitivity to very low neutrino fluxes. However, the relatively 
large spacing between sensors severely limits the sensitivity to lower energy muons (below 
100 GeV) which are the largest portion of the neutrino-induced muon spectra arising from 
annihilations in the Sun. To improve the sensitivity to these important lower-energy muons, 
a portion of the IceCube volume, referred to as DeepCore, has been outfitted with more 
densely packed sensors. 

The ability of IceCube/DeepCore to detect the flux from a given supersymmetric model 
depends on various factors such as the muon spectra, the neutrino backgrounds, and the 
method by which IceCube/DeepCore will analyze their results. However, to a very rough 
approximation, IceCube can detect muon fluxes on the order of 10 or 10^ /km^/yr above 
~100 GeV and DeepCore can detect muon fluxes on the order of 10^ or 10'^ /km^/yr above 
~10 GeV. In this paper, we primarily examine the latter case: total fluxes above 10 GeV. 
DeepCore can detect muons down ~10 GeV; however, the analysis threshold may turn out to 
be somewhat higher: Ref. jsof suggests an analysis threshold of ~35 GeV is reasonable, while 
Ref. [40] suggests 25-30 GeV is more likely (but possibly as low as 20 GeV) ^. In addition, 
the efficiency of detecting muons (or, equivalently, the effective area of the detector) falls 
with decreasing muon energy. Our results are not significantly affected by our choice of a 
10 GeV threshold as we use only an order-of- magnitude estimate of the IceCube/DeepCore 
sensitivity to the total muon flux above this energy. Our results are only affected if the flux 
is predominantly just above threshold {e.g. between 10 and 25 GeV), which is only expected 
to be the case for very light neutralinos (m^ <^ 100 GeV or mi/2 ^ 200 GeV) ^. 

We use the results of the WimpSim simulation jilj (as used within DarkSUSY j42[) to 
calculate the spectra of the neutrinos produced by the annihilations and the corresponding 
spectra of neutrino-induced muons in IceCube/DeepCore. More details of the neutralino 
capture/annihilation processes and the determination of the neutrino/muon fluxes may be 
found in Ref. 0]. 



III. REPRESENTATIVE STUDIES IN THE NUHMl 



A. Comparison with the CMSSM in the {mi/2,'mo) Plane 

The left panel of Figure |3] displays the neutrino-induced muon fluxes calculated in the 
{mi/2,mo) plane for the typical CMSSM scenario with tan/3 = 10, = and /i > 
introduced above. We recall that this plane has two narrow strips where the LSP density 
falls within the range allowed by WMAP and other measurements jol, which are coloured 
turquoise. One is the coannihilation strip close to the boundary of the forbidden charged- 



2 Triggering of two of IceCube/DeepCore's optical modules — which is possible for muons with energies 
as low as '-^lO GeV — is sufficient to detect a muon. However, two-module events suffer from extremely 
poor angular resolution. Muons that trigger three optical modules (which requires somewhat higher muon 
energies) yield much better track reconstruction and allow the analysis to be restricted to muons consistent 
with neutrinos coming from the direction of the Sun. 

The f GeV threshold was chosen and the bulk of the work done in this paper was performed prior to 
the availability of threshold estimates from Rcfs. 0,0- 

As this choice does not significantly affect our 
results, we have chosen to keep the current threshold in our analysis. 
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mi;2 (GeV) miy2 (GeV) 

FIG. 3: The (m'i/2)'^o) planes for tan 13 = 10 and Aq = in (left) the CMSSM and (right) the 
NUHMl for ji = 500 GeV. The solid (light blue) lines are contours of the neutrino-induced muon 
fluxes above 10 GeV in units of events /krn^ /yr. The shadings and other contours have the same 
meanings as in Figure 0. 

LSP region at low mo, and the other is the focus-point strip close to the EWSB boundary 
at large mo- 

As was discussed extensively in [3] and shown in the left panel of Figs. [T]&[2], at large mi/2 
it is usually not appropriate to assume that the LSP capture cross-section is dominated by 
spin- dependent interactions, nor that there is equilibrium between capture and annihilation. 
Accordingly, as already mentioned, neither assumption is made in this and subsequent plots. 
The solid (blue) lines are contours of the neutrino-induced muon fluxes above 10 GeV in 
units of events/km^ /yr. As seen in the left panel of Figure [3l the neutrino rate is potentially 
detectable in IceCube/DeepCore along a significant stretch of the focus-point strip, but 
apparently undetectable along the portion of the coannihilation strip that is compatible 
with the LEP Higgs constraint. As discussed in 0], the neutrino fluxes in the CMSSM are 
generally larger for tan /3 = 55 along the coannihilation strip (though still not observable in 
IceCube/DeepCore), and smaller along the focus-point strip than they are for the tan/3 = 10 
case shown. We also recall that when tan/3 = 55, there is also a third region compatible 
with the dark matter density constraint, namely a rapid- annihilation funnel at large mi/2, 
where the rate is far too small to be detectable by IceCube/DeepCore. 

The right panel of Figure [3] displays the neutrino-induced muon fluxes to be expected in 
a sample (mi/2, mo) plane in the NUHMl, again with tan /3 = 10 and Aq = 0, but now with 
fixed fi = 500 GeV. In this case, an excluded charged LSP region again appears at small mo 
though with a different shape from the CMSSM, and the EWSB boundary has moved to 
small mi/2 and mo. Here it is m\ that is driven negative rather than /i^ as in the CMSSM. 
There is a region favoured by — 2 that is almost excluded by the LEP Higgs constraint, 
and there is a coannihilation strip, as in the CMSSM, that follows the EWSB boundary. 
There is also a near-vertical extension of the coannihilation strip at mi/2 ~ 1000 GeV, that 
appears thanks to the freedom in the NUHMl of independently adjusting the /i parameter. 
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It occurs along a line where a particular relation between /i and mi/2 induces level crossing in 
the neutralino mass matrix as mi/2 is increased relative to /i, and thereby leads to increased 
Higgsino-gaugino mixing that brings the LSP density down into the WMAP range. To the 
right of this strip, the relic density falls below the WMAP range. This strip is the only part of 
this particular NUHMl plane where the neutrino- induced muon flux approaches detectability 
in IceCube/DeepCore, rising above 30 /km^/yr, thanks to the enhanced Higgsino-gaugino 
mixing. 

One of the key questions in our analysis will be the extent to which the freedom in the 
NUHMl, 2 to vary fi and/or provides this and other opportunities for IceCube/DeepCore 
detection of neutrinos that are absent in the CMSSM. 



B. NUHMl (^,mi/2) Planes 

We now illustrate further the behaviour of this IceCube/DeepCore-friendly WMAP strip 
with enhanced Higgsino-gaugino mixing, first in some representative {fi,mi/2) planes. Fig- 
ure m displays planes for tan (3 = 10 with (upper left) mo = 300 GeV and (upper right) 
mo = 500 GeV: both values of mo are in the range favoured by a frequentist analysis of the 



NUHMl parameter space [2l|. In each case we see EWSB boundaries at large and small 
mi/2 where m\ < 0, and regions with n < that are disfavoured by 6 — )■ sj. In the upper 
left panel, we also see charged LSP regions at large |yu| and mi/2- In each case, there is a pair 
of diagonal WMAP-compatible strips visible at ~ mi/2/2 where the Higgsino-gaugino 
mixing is enhanced, which are only weakly dependent on mo and are compatible with the 
LEP Higgs constraint for large enough mi/2- The neutrino-induced muon flux above 10 GeV 
is potentially detectable along essentially all of these diagonal WMAP-compatible strips in 
both panels, though decreasing as and mi/2 increase. These strips constitute exten- 
sions of the IceCube/DeepCore-friendly strip seen in the right panel of Figure [3l which has 
/i = 500 GeV and mi/2 ~ 1000 GeV, to different values of these NUHMl parameters. The 
red dot-dashed curve is the contour for ruh = 114 GeV and one should preferably lie above 
this curve, though one should recall that there is a 1.5 GeV uncertainty in the theoretical 
calculation of ruh- Nevertheless, fluxes above the Higgs limit still reach above 500 /km^/yr. 
For mo = 300 GeV, the g — 2 constraint would prefer lower values of mi/2 in potential conflict 
with the Higgs bound - though the muon fluxes are quite large where both constraints are 
satisfied. At mo = 500 GeV, the g — 2 is not satisfied within 2a anywhere on the plot. We 
note that there is another WMAP strip slightly above the region with no EWSB. This strip 
corresponds to the heavy Higgs funnel where 2m^ ^ m^ and would not be present in the 
CMSSM at tan /3 = 10. However, neutrino-induced muon fluxes are too small along these 
strips to be observed in IceCube/DeepCore. 

Similar diagonal strips with enhanced Higgsino-gaugino mixing are seen in the lower left 
panel of Figure H] for tan /3 = 20 and mo = 500 GeV, though the muon fluxes are generally 
smaller. However, the diagonal strip for /x > has a part where — 2 lies within the 
favoured range and the muon flux is > 300 /km^/yr. As in the previous panels of Figure HJ 
the WMAP-compatible strips corresponding to the funnel and close to the EWSB boundary 
have muon fluxes which are unobservably low. In the lower right panel of Figure H] for 
tan/3 = 55 and mo = 500 GeV we see that, whilst the EWSB boundary is located similarly 
to the previous panels, the charged LSP region has advanced considerably, as has the region 
excluded by 6 — )■ 57. At this value of tan/3, /x < is not consistent with RGE running. 
We also note that the diagonal strip with enhanced Higgsino-gaugino mixing has merged 
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FIG. 4: The (/^, 7711/2) planes in the NUHMl for ^0 = and (upper) tan/3 = 10, (left) rriQ = 
300 GeV and (right) ruo = 500 GeV, (lower) mo = 500 GeV and (left) tan/3 = 20 and (right) 
tan/3 = 55. The shadings and contours have the same meanings as in Figurel^i 



with the rapid-annihilation funnel, with both sides of the funnel now easily discernible. The 
muon fluxes where 6 — )■ S7 is acceptable are <50 /km^/yr and large enough to be detectable 
by IceCube/DeepCore in the lower part of the transition strip, for values of mi/2 above the 
bounds imposed by the LEP Higgs constraints and well with the la bounds from g — 2 
(denoted by the dashed black curves). The entire region with mi/2 below 700-800 GeV 
yields values oi — 2 in the range favoured by experiment. In contrast to the cases with 
lower tan /3, the funnel region now also has neutrino fluxes that are sufficiently large to be 
observable ^. 



We note that equilibrium is established along the transition strip in all cases, but not along the funnel 
except at very large tan/3. We also note that spin-dependent scattering is dominant at low |//| for 
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C. NUHMl (/i,mo) Planes 

In order to explore further the IceCube/DeepCore-friendly region, in Figure [5] we display 
NUHMl planes for tan/3 = 10 with (upper left) mi/2 = 300 GeV and (upper right) mi/2 = 
500 GeV (again, both values of mi/2 are in the range favoured by a frequentist analysis of 
the NUHMl parameter space [sH), and for mi/2 = 500 GeV with (lower left) tan /3 = 20 and 
(lower right) tan /3 = 55. We see in both the upper panels EWSB boundaries at large and 
small mo, charged LSP regions at smaller and mg, and regions excluded by 6 — )■ 57 when 
/i < 0. When mi/2 = 500 GeV, there is a small region (shaded black) between the stau LSP 
region and the EWSB boundary where the LSP is a right-handed selectron (or smuon). As 
in the right panel of Figure [3l we see in both of the upper panels WMAP-compatible strips 
at roughly fixed values of related to the values of mi/2, along which the muon flux above 
10 GeV is potentially detectable: above 500 events/km^/yr for (upper left) mi/2 = 300 GeV 
and /i ~ 200 GeV, and above 100 events/km^/yr for (upper right) mi/2 = 500 GeV and 
/i ~ 300 GeV. As already commented, along these strips the fluxes are enhanced by Higgsino- 
gaugino mixing. We note that only part of the strip for mi/2 = 300 GeV is compatible with 
the LEP Higgs constraint, whereas all the mi/2 = 500 GeV strip is compatible. There are no 
parts of the IceCube/DeepCore- friendly strips in regions favoured by (7^ — 2. In both panels, 
there are WMAP-compatible extensions of these strips to larger /i at uiq ~ 100 GeV, due 
to coannihilation, which segue into funnel strips close to the EWSB boundary. However, 
these do not yield neutrino fluxes interesting for IceCube/DeepCore. In both upper panels, 
equilibrium is established and spin-dependent scattering is dominant along the transition 
strip, whilst the opposite is true along the funnel. 

Turning now to the lower panels in Figure we see that as tan /3 increases with fixed 
mi/2 = 500 GeV the EWSB boundary moves away to larger the charged LSP region rises 
to larger mo and (for tan /3 = 55) the 6 — )■ 57 exclusion extends to /i > as well as the visible 
parts of the half-plane with /i < 0. In the lower left panel for mi/2 = 500 GeV and tan /3 = 20 
we see near- vertical WMAP-compatible strips with ~ 300 GeV where Higgsino-gaugino 
mixing is enhanced and the neutrino flux is favourable for IceCube/DeepCore, and part 
of the strip for > is also favoured by — 2. This interesting region then bends 
into a near-horizontal coannihilation strip, resembling those in the upper panels, where the 
neutrino-induced muon flux is mostly unfavourable. As |/i| increases at low mo, the LSP 
changes from a mostly right-handed stau to a right-handed e/jx. The muon flux is also 
small in the funnel that follows the stau LSP boundary. Turning finally to the lower right 
panel in Figure O we clearly see a two-sided diagonal funnel where the rapid annihilation 
via direct-channel heavy Higgs poles brings the relic density into the WMAP-compatible 
range, albeit with a relatively low muon flux. We also see that the strip with enhanced 
Higgsino-gaugino mixing is less vertical and with a lower flux than previously, though still 
IceCube/DeepCore-friendly. In this case, the g — 2 constraint is satisfied over much of the 
plane. 



tan/3 = 10 and 20, but subdominant at larger for all values of tan/3. 
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FIG. 5: The (/i,mo) planes in the NUHMl for Aq 
300 GeV and (right) mx/2 = 500 GeV, (lower) 
tan/3 = 55. In the black regions the LSP is the right-handed selectron or smuon, and the other 
shadings and contours have the same meanings as in Figure\^ 



and (upper) tan/3 = 10, (left) mi/2 = 
500 GeV and (left) tan 13 = 20 and (right) 



D. NUHMl {mA,mi/2) Planes 



Figure [6] displays some representative {mA,iTii/2) planes in the NUHMl for Aq = and 
tan/3 = 10 with (upper left) mo = 300 GeV and (upper right) mo = 500 GeV. In both cases, 
we see that the EWSB requirement excludes a triangular region at large and small mi/2 
where /x^ is driven negative, whereas 6 — )■ 57 excludes a band at small m^. There is a region 
favoured by g^~2 in the left plane that is again almost excluded by the LEP Higgs constraint. 
In each case, a WMAP-compatible strip runs parallel to the EWSB boundary, extending to 
small TTiA at mi/2 ~ 125 GeV. There is also a diagonal funnel where LSPs annihilate rapidly 
though direct-channel Higgs poles, because ~ mA/2, and there are WMAP strips on 
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FIG. 6: The (mA,rni/2) planes in the NUHMl for Aq = and (upper row) tan/3 = 10 and (left) 
mo = 300 GeV, (right) rtiQ = 500 GeV, (lower row) for mo = 500 GeV and (left) tan/3 = 20 and 
(right) tan/3 = 55. The shadings and contours have the same meanings as in Figure\M 

either side of this funnel. The neutrino flux is very unfavourable along the funnel, since the 
neutrino-induced muon flux above 10 GeV is <^ 10 event s/km^/yr. However, a signal may 
be observable in IceCube/DeepCore along the EWSB boundary, where the muon flux above 
10 GeV lies in the range 10 to 100 events/km^ /yr. As we have seen before, equilibrium is 
established and spin-dependent scattering is dominant along the transition strip near the 
EWSB boundary, in contrast to the funnel regions where equilibrium is not established and 
spin-independent scattering is dominant. 

Turning now to the lower row of plots in Figure [6] for (left) tan /3 = 20 and (right) 
tan /3 = 55, both with mo = 500 GeV, we see that the EWSB boundary is broadly similar, 
whereas the region forbidden by 6 — )■ 57 is much reduced. The WMAP-compatible rapid- 
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annihilation funnel is clearly visible for tan /3 = 20, but folds into the EWSB boundary 
strip for tan /3 = 55. As in the previous cases, the neutrino flux is unobservably small along 
the rapid-annihilation funnel, but may be observable along the EWSB boundary strip for 
tan (3 = 20. We note that a portion of the IceCube/DeepCore-friendly region for tan /3 = 20 
is favoured by — For tan /3 = 55, muon fluxes are > 10 /km^/yr for mi/2 500 GeV along 
the transition/funnel and compatible with g — 2. We note that this case is also interesting 
because in much of the allowable plane capture is dominated by spin-independent scattering. 

E. NUHMl {mA,mo) Planes 

Figure [7] shows some sample NUHMl (m^, mg) planes with fixed tan /3 that exhibit rapid- 
annihilation funnels. In these cases, the funnels appear as essentially vertical double strips on 
either side of the line where = 2m^. In each case, we also note the presences of WMAP- 
compatible strips close to the EWSB boundary where fi^ = 0. The two are attached by 
coannihilation strips at mo ~ 100 GeV. When mi/2 = 500 GeV, coannihilation to the left of 
the funnel at low ttla is dominated by selectrons/smuons which are the LSPs in the lower left 
corners. In the upper left panel for tan /3 = 10 and mi/2 = 300 GeV, we also see a region at 
low rriQ that is favoured by — but this is in a region disfavoured by the LEP Higgs limit. 
In this plane, the only WMAP-compatible region allowed by the other constraints is up the 
funnel at large mo, where the neutrino flux in unobservably low. In the upper right panel for 
^1/2 = 500 GeV, the Higgs constraint is irrelevant, but the neutrino flux is still very low in 
all the WMAP-compatible region except along the strip close to the EWSB boundary. The 
same is true in the lower left plot for tan/3 = 20 and mi/2 = 500 GeV, where we note that 
a part of this strip is inside the region favoured hy g^ — 2. However, as shown in the lower 
right plot, when tan/3 is increased to 55 with the same value of mi/2, the EWSB boundary 
moves so close to the rapid-annihilation funnel that there is no IceCube/DeepCore-friendly 
region. 

F. NUHMl Summary 

We have found in the above analysis that interesting neutrino fluxes may arise from 
enhanced Higgsino-gaugino mixing. This possibility arose in the CMSSM in the focus-point 
region along the EWSB boundary, but may arise in regions of the NUHMl parameter space 
that are far from this boundary, corresponding in general to larger values of mi/2- The 
EWSB boundary region may also yield interesting fluxes when /i^ — )■ at the boundary. 
One of the other mechanisms present in the NUHMl for bring the relic LSP density into the 
WMAP-compatible range is rapid annihilation through direct-channel Higgs poles. However, 
we have found that this region generally yields low neutrino fluxes. 

IV. REPRESENTATIVE STUDIES IN THE NUHM2 

As discussed earlier, in the NUHM2 we are free to adjust both soft Higgs masses inde- 
pendently and as a consequence we can study models for which both /i and niA are free 
parameters [iS]. We flrst consider some {fi,mA) planes for Aq = 0,mi/2 = 300 GeV and 
mo = 100 GeV as shown in Figure [H In each case, we see a WMAP-compatible rapid- 
annihilation funnel at ttia ~ 250 GeV as the neutralino mass is roughly 0.43 mi/2 and 
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FIG. 7: The (m^,mo) planes in the NUHMl for Aq = Q and (upper row) tan/3 = 10 and (left) 
mi/2 = 300 GeV, (right) mi/2 = 500 GeV, and (lower row) for = 500 GeV and (left) 

tan/3 = 20, (right) tan/3 = 55. The shadings and contours have the same meanings as in Figure\^ 



therefore roughly half the pseudoscalar mass across the plane. In addition there is an arc at 
relatively large /i and rriA due to the coannihilations between the neutralino and sneutrinos. 
The dark blue shaded regions have a sneutrino LSP and are excluded 4^. In the left panel 
for tan/3 = 10, we also see a near- vertical strip at low where the LSP has a substantial 
Higgsino component, which extends downwards in a stau coannihilation strip almost parallel 
to the charged dark LSP boundary. The > half of the plane is consistent with the g — 2 
constraint. As we have seen before, equilibrium is in this case well established along the ver- 
tical transition strip, partially established along the horizontal funnel, and not at all along 
the arcing coannihilation strip. While scattering is heavily dominated by the spin-dependent 
cross-section along the transition strip, spin-independent scattering is significant along the 



16 




FIG. 8: The {n,mA) planes in the NUHM2 for Aq = 0,mi/2 = 300 GeV and mo = 100 GeV with 
(left) tan/3 = 10 and (right) tan/3 = 55. In the blue shaded regions the LSP is a sneutrino, and 
the other shadings and contours have the same meanings as in Figure\^ 

other two strips. In the right panel for tan /3 = 55, the charged LSP boundary is present at 
lower /i and more prevalent at large itia- Both the stau and sneutino LSP regions are again 
parallelled by WMAP-compatible coannihilation strips. The supersymmetric contribution 
to 5^ — 2 in the right panel exceeds the 2a bound. 

In the left plane, the only region with an IceCube/DeepCore-friendly neutrino flux is the 
vertical strip with a substantial Higgsino content, analogous to those seen previously in our 
NUHMl analysis. Here fluxes are in excess of 500 /km^/yr. On the other hand, in the right 
panel where this strip has disappeared, the only IceCube/DeepCore-friendly region is at 
small /i and niA-, below the rapid-annihilation funnel. We further note that, in both planes, 
the Higgs mass falls below 114 GeV. 

Similar trends can be seen in Figure [9] for Aq = 0,mi/2 = 500 GeV and mo = 300 GeV. 
In the left panel for tan/3 = 10, we see a near-vertical strip where the LSP has an enhanced 
Higgsino component, split in two by a near- horizontal rapid-annihilation funnel. This funnel 
is also present in the right panel for tan/3 = 55, and splits a diagonal coannihilation strip 
close to the charged LSP boundary. Again, as in the previous case for smaller mi/2 and 
mo, the only IceCube/DeepCore-friendly regions are in the Higgsino-like strip at small n 
for tan /3 = 10, and for very small fi and rriA for tan/3 = 55. Spin-independent scattering 
becomes dominant for |yu|>500 GeV when tan /3 = 10 and is dominant everywhere in the 
displayed plane for tan/3 = 55. In this case, the Higgs mass is above 114 GeV everywhere 
above the red dot dashed curved found at low tua- For tan/3 = 10, the supersymmetric 
contribution to (7 — 2 is too small across the plane, whereas for tan /3 = 55 it lies within the 
experimental range. 

Finally, we present some {1711,1712) planes in Figs. [TO] and [TT| for mi/2 = 300 GeV, mo = 
100 GeV and mi/2 = 500 GeV, mo = 300 GeV, respectively. In each case, we assume Aq = 
and the left panel is for tan /3 = 10 and the right panel for tan/3 = 55. The ranges of mi 
and m2 are symmetric G (—1000, 1000) GeV, except for the right panel of Figure [TU| where 
asymmetric ranges (mi, m2) G (—2000, 0) GeV are chosen so as to display better the regions 
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FIG. 9: The {fi, mo) planes in the NUHM2 for Aq = 0, mi/2 = 500 GeV and mo = 300 GeV with 
(left) tan/3 = 10 and (right) taii/3 = 55. The shadings and contours have the same meanings as in 
Figure 
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FIG. 10: The (mi, 7712) planes in the NUHM2 for mi/2 = 300 GeV and mo = 100 GeV with (left) 
tan 13 = 10 and (right) tan^S = 55. The signs shown for mt actually refer to the sign of mf as they 
are run in the RGBs. The shadings and contours have the same meanings as in FigurelEi 

not excluded by the EWSB constraint (pink) and the sneutrino LSP regions (blue) or the 
charged-LSP constraint (brown). The signs shown for rrii actually refer to the sign of mf as 
they are run in the RGEs. As these are GUT scale parameters, negative values may indicate 
a cosmological issue with broken symmetric vacua. For a recent discussion of this issue in 
the NUHM, see 0. 

In the left panel of Figure [TUl for tan /3 = 10, we see two narrow WMAP strips following 
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the top left and bottom right boundaries of the allowed lozenge of the (mi, 7712) plane. Of 
these, the former has an IceCube/DeepCore- friendly neutrino-induced muon flux over most 
of its length, rising to > 300 events/km^/yr for m2 ~ 400 GeV. Note that the nearly hor- 
izontal part of the strip corresponds again to the transition strip where a there significant 
Higgsino contribution to the neutralino composition. The more vertical part of the strip 
corresponds to the rapid annihilation funnel. On the other hand, the bottom right WMAP 
strip corresponding to neutralino-sneutrino coannihilation generally has an unobservably 
small neutrino flux, except for positive values of m2, where the muon flux may reach 30 
events/km^ /yr. As usual, the larger neutrino fluxes are reached when the Higgsino compo- 
nent of the LSP is enhanced. In this plane the g ~2 constraint is satisfied, though the Higgs 
mass is low. 

In the right panel of Figure [TOl for tan/3 = 55, the allowed region of the (mi, 777,2) plane 
has receded to more negative values of mi and m2, and the two narrow WMAP strips are 
squeezed closer together. Additionally, we note that the top left strip has bifurcated along 
the two sides of a rapid-annihilation funnel, whereas the lower right strip follows either the 
stau or sneutrino coannihilation boundaries. As in the left panel of Figure [101 the neutrino 
fluxes are generally more favourable along the top left strip, though they also become more 
IceCube/DeepCore-friendly towards the upper end of the other strip, near their junction. 
Neither g — 2 or the Higgs mass constraints are satisfied in this plane. 

In the left panel of Figure [TT] for tan/3 = 10, we see a single narrow WMAP strip following 
the top boundary of the allowed region of the (mi, m2) plane corresponding to the transition 
strip. As in the corresponding panel of Figure [TUl the neutrino- induced muon flux is largest, 
exceeding 100 events/km^/yr and hence quite IceCube/DeepCore-friendly, in the top part of 
the strip close to the EWSB boundary where m2 ~ 800 GeV and the LSP has an enhanced 
Higgsino component. The left boundary corresponds to the the funnel region where there 
are two strips, one of which continues on to the vertical transitions strip. The upper part 
of the funnel still has observable fluxes, but these quickly drop as m2 is decreased. Here, 
and in the right panel as well, equilibrium is established along all WMAP strips, however 
spin- dependent scattering is sub-dominant almost everywhere in the plane. 

In the right panel of Figure [TT] for tan /3 = 55, the previous single narrow WMAP strip 
has again bifurcated along the two sides of a rapid- annihilation funnel, and has moved away 
from the boundaries of the allowed region of the (mi,m2) plane. Once again, the neutrino 
flux is largest for m2 > 0, though generally lower than in the left panel of Figure [TT] (where 
tan /3 is smaller) or in the right panel of Figure [TO] (where mi/2 and tuq are smaller), and only 
barely IceCube/DeepCore-friendly. In this plane, both the g — 2 and Higgs mass constraints 
are satisfied. 

V. SUMMARY 

In our previous analysis of the CMSSM, we found that the flux of high-energy neutrinos 
from LSP annihilations inside the Sun was likely to be observable along the focus-point 
WMAP strip, where the Higgsino component of the LSP is enhanced, and at the tip of the 
coannihilation WMAP strip where the LSP is relatively light. On the other hand, there were 
significant portions of the WMAP-compatible strips in parameter space, particularly along 
the coannihilation strip and heavy-Higgs rapid-annihilation funnels, where the the neutrino 
flux was not IceCube/DeepCore-friendly. 

We find some similar features in our analyses of the NUHMl and NUHM2. Specifically, 
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FIG. 11: The (mi, 7712) planes in the NUHM2 for mi/2 = 500 GeV and mo = 300 GeV with (left) 
taii/3 = 10 and (right) tan/3 = 55. The shadings and contours have the same meanings as in 
Figure O 

there are significant portions of tlie WMAP strips where the high-energy solar neutrino flux 
is unhkely to be observable with IceCube/DeepCore. In these models, IceCube/DeepCore- 
friendly fluxes are often found in regions where the LSP has an enhanced Higgsino compo- 
nent, and this occurs under circumstances that cannot be realized in the CMSSM. Specifi- 
cally, it may occur for larger LSP masses than along the focus-point strip of the CMSSM: 
see, for example, the right panel of Figure [3l the first three panels of Figure HI and the first 
three panels of Figure [51 

We conclude, therefore, that IceCube/DeepCore has interesting prospects for probing 
aspects of the NUHMl and NUHM2 parameter spaces. However, it seems clear that a 
more complete exploration of these models, capable of measuring a high-energy neutrino 
flux when the LSP is relatively heavy and/or does not have a large Higgsino component, 
would require a subsequent generation of experiment. On the other hand, one may hope 
that forthcoming LHC results and/or direct searches for LSP scattering could provide more 
encouraging indications on the prospects for searches for supersymmetric dark matter via 
annihilation into high-energy solar neutrinos. 
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